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Abstract 

We discuss the application of a variant of the method of simplest 
equation for obtaining exact traveling wave solutions of a class of non¬ 
linear partial differential equations containing polynomial nonlinear¬ 
ities. As simplest equation we use differential equation for a special 
function that contains as particular cases trigonometric and hyperbolic 
functions as well as the elliptic function of Weierstrass and Jacobi. We 
show that for this case the studied class of nonlinear partial differen¬ 
tial equations can be reduced to a system of two equations containing 
polynomials of the unknown functions. This system may be further re¬ 
duced to a system of nonlinear algebraic equations for the parameters 
of the solved equation and parameters of the solution. Any nontriv¬ 
ial solution of the last system leads to a traveling wave solution of 
the solved nonlinear partial differential equation. The methodology 
is illustrated by obtaining solitary wave solutions for the generalized 
Korteweg-deVries equation and by obtaining solutions of the higher 
order Korteweg-de Vries equation. 
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1 Introduction 


Nonlinearity is an essential feature of many systems in Nature and society 

ra-ra 

Traveling wave solutions of nonlinear partial differential equations are 
studied much in the last decades [ITj-|2Ij as they occur in many natural sys¬ 
tems [22]-[25] and because of existence of various methods for obtaining such 
solutions [26] - [30]. Below we shall consider the method of simplest equation 
for obtaining exact analytical solutions of nonlinear partial differential equa¬ 
tions [31]- [M] and especially its version called modified method of simplest 
equation [35] - [37]. Method of simplest equation is based on a procedure 
analogous to the first step of the test for the Painleve property [31], [55] . 
In the version of the method called modified method of the simplest equa¬ 
tion [59]. m this procedure is substituted by the concept for the balance 
equation. Modified method of simplest equation has its roots back in the 
history (for an example see [H]-|ll]). Method of simplest equation has been 
successfully applied for obtaining exact traveling wave solutions of numerous 
nonlinear PDEs such as versions of generalized Kuramoto - Sivashinsky equa¬ 
tion, reaction - diffusion equation, reaction - telegraph equation, generalized 
Swift - ffohenberg equation and generalized Rayleigh equation, generalized 
Fisher equation, generalized ffuxlcy equation , generalized Degasperis - Pro- 
cesi equation and b-equation, extended Korteweg-de Vries equation , etc. 

051 - 0 - 

A short summary of the method of simplest equation is as follows. First 
of all by means of an appropriate ansatz (for an example the traveling-wave 
ansatz) the solved of nonlinear partial differential equation is reduced to a 
nonlinear ordinary differential equation 

P(u,u f ,u K ,...) = 0 (1.1) 

Then the finite-series solution 

«K) = jt !V[»(«)]“ (1-2) 

fl= — u 

is substituted in (EUi. p M are coefficients and g(£) is solution of simpler 
ordinary differential equation called simplest equation. Let the result of this 
substitution be a polynomial of <?(£)• Eq. (jl .2]) is a solution of Eq. 01.ip if all 
coefficients of the obtained polynomial of g(£) are equal to 0. This condition 
leads to a system of nonlinear algebraic equations. Each solution of the 
last system leads to a solution of the studied nonlinear partial differential 
equation. 
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In this article we consider a large class of (l+l)-dimensional nonlinear 
partial differential equations that are constricted by polynomials of the un¬ 
known function and its derivatives. As simplest equation we shall use equa¬ 
tion of the kind 



The text below is organized as follows. In Sect. 2 we introduce the class of 
studied nonlinear partial differential equations and the used class of simplest 
equations and their solutions. Then we show that any of the nonlinear par¬ 
tial differential equations of the discussed class can be reduced to a system 
of two equations containing polynomials of the unknown function. These 
polynomials can be obtained on the basis of addition and multiplication of 
some basic polynomials connected to the derivatives of the solved nonlinear 
partial differential equation. In section 3 we calculate some of the most used 
basic polynomials. In Sect. 4 the methodology is illustrated by application 
for obtaining solitary wave solutions of 

• generalized Korteweg-deVries equation; 

• higher order Korteweg-deVries equation. 

Several concluding remarks are summarized in Sect. 5. 


2 Formulation of the method 


2.1 Proof of the basic theorem 


Let us consider a nonlinear PDE with nonlinearities that are polynomials of 
the unknown function h(x,t) and its derivatives. We search solution of the 
kind 

h(x,t) — /i(£); £, = H'X + vt (2.1) 

where /i and v are parameters. The basis of our search will be a solution 
g(£) of a certain simplest equation. Hence 


h = /[»(£)] 


( 2 . 2 ) 


h from Eq. (12. 2 H is a composite function. For the n-th derivative of h we have 
the Faa di Bruno formula [52] 


h(n) 


5Z/(fc) 


k =1 


p(fc,n) 


n\ 


A i 

%) 


n (A,!)(i! 




(2.3) 


where 
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7 _ d n h . 

'HO dx n ’ 


f = dk f ■ 

J (fc) dgk ) 


_ d l g . 

^(*) — da:*’ 


p(n, k) = {Ai, A 2 , • • •, A n }: set of numbers such that 

n n 

y \ = k-,y i A* 


= n. 


(2.4) 


i= 1 




Further we shall concentrate on fn^ and gg\. 

Let us now assume that / is a polynomial of g. Then 

/ = X>9 r 

r =0 


(2.5) 


Let us consider the derivative fgg. If k > r this derivative is 0. The derivative 
is non-zero if k < r. We shall use the function (~) r k with the following 
definition 

^ f 0 r < k _ „ s 

Qrk = < ( 2 . 6 ) 

II r > k 

Then the derivative fgg can be written as 

<? 


f(k) ^ ^ 


T ! 


r=0 


(r — k)\ 


b r g 


r—k 


(2.7) 


The derivative gg\ is connected to the simplest equation. In general we can 
use the following simplest equation with polynomial nonlinearity 


9{k)= (4f) =y a :9- 

^ ' 3=0 


3 


( 2 . 8 ) 


The solution of this equation defines certain special function I4 0)aii ... iam (£; k, /, m) 
where 


• k: order of derivative of g\ 

• l: degree of derivative in the defining ODE; 

• m: highest degree of the polynomial of g in the defining ODE. 
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This special function has very interesting properties as its particular cases 
are the trigonometric, hyperbolic, elliptic functions of Jacobi, etc. Below we 
shall use the function I4 0)aii ... iam (£; 1,2, m) which is solution of the simplest 
equation 



We let m undetermined for now. 

If g 2 ^ is given by Eq. fl2.9p then what is the relationship for the derivative 
fi'h ) 7 

Lemma. If g 2 L) i s given by Eq. 1 2.91) then the following relationship holds 
for the derivative g^ : 


9(f) = [A(g)} Xz O itXi (g, g(i))fi\ t {9(i)) 

where 


• \i is a non-negative integer; 

• Ai(g): polynomial of g; 


Oi,Xi{g,g(i)) 


i — even , Aj — even 
i — even, A i — odd 

< a j9 j Y Xi ^ i — odd, Xi — even 

3=0 

m 

(XI i — odd, Xi — odd 

. 3=0 



9 ( 1 ) 


m 


where [Aj/2] denotes the integer part o/A*/2; 


1 )) 


1 = i — even 

g; i) i — odd 


Proof. It is easy to show (by direct differentiation and by induction) that if 
g 2 X] is given by Eq. (12.91) then 


9{2n) — A 2 n(fiO; 9{2n+l) ~ A 2n+ i (<7)<7(1) (2.10) 

where A 2 n (g) and A 2 n+ 1 (g) are polynomials of g. 
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Let then 


1 = -o 


. 9 ( 1 ) * - even 

9(i) * - odd 


( 2 . 11 ) 


^i(9(l)) 

Hence Eq. (j2.10p can be written as 

9(i) = A(g)tti(g { i)) (2.12) 

Thus 

9( A .) = [■ 4 i(9)] A ‘P.(j(i))] A ‘ (2.13) 

Now from the definition of - Eg. (12. lip we see that is equal to 1 if i is 

even. is equal to gA if i is odd. But from Eq. (12.101) for the case of odd 

A,: 


n \ _ 2([Aj/2])+l _ 

%) “ %) 


= 9(i)(j ( 2 i,) IA ‘ /21 (2.14) 


In addition from Eq. (12.101) it follows that is a polynomial of g (this also 
can be seen from the definition of the simplest equation (12.9K ). Thus we can 
define 


1 

i 

3(1) 


i — even, A; — even 
i — even, Aj — odd 


Oix(g>9( i)) = < 


(E a jg j )^ i — odd, Xi 


even 


(2.15) 


l-o 

m 


(Z) a jg j y Xi ^ i — odd, A i — odd 
j =o 


(2.16) 


and then 

[ fi i(9(i))] Ai = Oi^(g,g {1) )nxi (9(i)) 

Substitution of Eq. (l2.16p in Eq. (l2.13p leads us to the relationship we want 
to prove 

4 = [-4.(9)] A -O iA ( J , 9a) )ll A ,( J(1) ) (2.17) 

□ 


Now we are in position to prove 

Theorem. If g^ is given by Eq. Id. ,91) and f is a polynomial of g given by 
Eq. / Id. 51) then for h[f(g)] the following relationship holds 

h(n) = Kn(q,m)(g) + g{i)Z n {q, m)(g) 

where K n (q,m)(g ) and Z n (q,m)(g ) are polynomials of the function c/(£). 


6 


















Proof. The substitution of Eqs. (12.17j) and (12.7j) in Eq. (l2.3p leads to the 
following relationship for 


n q 


k =1 r =0 ' ' p(n,k) i =1 


(r — /c)! ^ 'J-l (Aj!)(i!) A 

V ’ p(n,k) i= 1 V WV y 

Now let us show that the relationship (12.181) can be written as 


(2.18) 


ft(„) = K n (q,m)(g) + gii)Z n (q,m)(g) (2.19) 

where K n (q,m)(g ) and Z n (q,m)(g ) are polynomials of the function g(£). I 11 
oder to show this we rewrite Eq. (12.181) as 


n q ; 

fc=l r =0 ' 


fc)! 


&r<7 r “* ^ n\ 


p(n,k) 


fritter 

11 (A,!)(i!)** 


( 2 . 20 ) 


2—1 


and consider P = n <^1,A* , P(i))(^(i)) ■ T llis product is equal of polyno- 
2=1 

rnial of g multiplied by g"^ where cr is an integer. There are two possibilities 
for a: 


• a: even . Then (gj^) 0 '/ 2 is a polynomial of g according to Eq. (12.91) 

• a: odd. Then g'E is equal to <7(i)(<7(i))^ 2 and according to Eq. (12.91) 
this is equal to g^ multiplied by a polynomial of g. 

Because of all above there will be two kinds of terms in h( n y. terms that are 
polynomials of g and terms that contain gm multiplied by a polynomial of 
g. Collecting the two kinds of terms we arrive at Eq. (|2.19p . □ 

Let us note that for some values of n one of the polynomials K n (q , m) or 
Z n (q,m ) can be equal to 0. 


2.2 Formulation of the method 

On the basis of all above the modified method of simplest equation based on 
simplest equation (12. 9 p can be formulated as follows: 

1. We consider a nonlinear PDE E*(u(x, t),u x (x, t),u t (x, t),u xt (x, t ),...) = 
0 where E* is a polynomial of u(x, t) and its derivatives. We search for 
solutions of this equation based on the ansatz £ = fix + vt where a and 
fd are parameters. 
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2. The ansatz £ = /j,x + vt reduces the nonlinear PDE to the ODE 
E(h, /i(i),...) = 0 where E is a polynomial of h and its derivatives. 

3. We assume h and g are given by the Eqs. (12.51) and (I2.9[) . Substitution 
of the relationships in the equation E — 0 reduces any derivative of 
this equation to a term of the kind (j2. 19j) . 

4. As the terms in E — 0 are polynomials of h and its derivatives then 
the equation reduces to a polynomial containing g , g^, g^, ...: E = 

^o*(^)+^ / r(5 , )5 , (i) + ^ 2 *( 5 , ) 5 , p) +_Eq. (12.90 reduces the higher degrees 

(n > 1 ) of ( 7 ( 1 ) to a polynomial of g (for even n) or to a polynomial of 
g multiplied by g (!) (for odd n). Thus equation E = 0 is reduces to 

E = W Q (g) + Wi(g)g(i) = 0 (2.21) 

where 1 Toy (( 7 ) are polynomials of g. 

5. In order to obtain a nontrivial solution of Eq. (l2.21|) we have to balance 
the highest degree of the polynomial Wq (i.e., to ensure that there are 
at least two terms that contain the highest degree of Wo). The same 
has to be made for the polynomial Wi. As a result we obtain one or two 
relationships among the parameters of the equation and parameters of 
the solution. These equations are called balance equations. Balance 
equations fix parameters q and m. 

6 . Further we set to 0 all coefficients of the polynomials W 0 (g) and W\{g). 
The result is a system of nonlinear algebraic equations that contains 
the parameters of the equation, parameters of the solution ( 12 .5|) and 
parameters of the simplest equation (12.91) . 

7. Any nontrivial solution of the above system of algebraic equations (if it 
exists) leads to a traveling wave solution of the nonlinear PDE E* = 0. 


Calculation of some of polynomials K n and 
Z n from Eq.( l2.19F ) 


The derivatives /rp), /q 2 ), h( 3 ), /i( 4 ), /q 5 ), /?.( 6 ) and h ( 7 ) are much used in the 
model nonlinear partial differential equations. Below we shall calculate the 
polynomials K n and Z n connected to these derivatives. 

Let us mention first that the polynomials K n+ 1 and Z n+ 1 are connected 
to the polynomials K n and Z n . This relationship can be obtained on the 
















basis of the relationship 


(3.1) 


/ d , 

V+l) = ^(n) 

Substitution of Eg. (12. 19(1 in Eg. (13.Ill leads to 


dZ n 2 
Zn9m + "di -9 '" 


d(n+ 1) 

Taking in account Eg. (l2.9p we obtain 
Kn +1 

^ 71+1 


+ 


dKr, 

dg 


-9(i) 


Y m m 

■ 7_i . X - '' j 

+ -^-E‘w’ 

i=o y i=o 

dKji 


dg 


(3.2) 


(3.3) 


Hence we need to calculate only K] and Z\ and then we can obtain K n and 
Z n , n = 2, 3,... by the recurrence eguations (13.3H . 

We can write 


K 0 

Z 0 


J2 b -s r 

r =0 


0 


From Eg. (12.1811 we obtain 


Then 


Ki — 0; Z\ — rb r g r ~ l 


r=0 


q m 


= EE 

r=0 j=0 L 

•^2 = 0 . 


jr + r(r — 1) 


b r Q j+r ~ 2 


CLjU r Q' 


K, 

Z 3 


0 ; 

q m 

EE 

r=0 j =0 L 


—jr + r(r — 1) 


(j + r - 2 )a j b r g ,+r 3 


(3.4) 


(3.5) 


(3.6) 


(3.7) 
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q m m r- 


* = EEE 

r=0 j =0 u =0 

Z 4 = 0. 


Ij r _|_ r(r - l)\ (j + r - 2) Qtt + j + r - 3 


ajb r a u g-’ +r+u ~ 4 

(3.8) 


K b 

e 5 


= 0 ; 


q m m 


EEE 

r=0 j=0 w=0 L 

r + u — 4 )cijb r a u g :j+r+u ~ 5 


Ij r + r ( r _ i)^ (j + r ~ 2 ) + j + r - 3 


(j + 

(3.9) 


IU = 


q m m m y 

EEEE 

r=0 j=0 w=0 v=0 


lj r + r (r - 1)^ (j + r - 2) + j + r - 3 ) (j + 


Z fl = 


r + u — 4) 

0. 


(^v + j + r + u — 5^ ajb r a u a v g j+r+u+v 6 


(3.10) 


K r 

Z 7 


0; 

q m m m y 

EEEE 

r=0 j= 0 w=0 


jr + r (r - 1) (j + r - 2) ( -u + j + r - 3 (j + 


r + u — 4) 


-v + j + r + w- 5j (j + r + u + v - 6)ajb r a u a v gj +r+u+v 7 

(3.11) 


etc. 

For the practical application of the modified method of simplest equation 
we need to calculate the maximum grade of polynomials in h( n y As we have 
seen above in the text /q n ) consists of two kinds of terms: polynomial of 
g plus another polynomial of g multiplied by g^y We note that the above 
maximum grades have to be non-negative. Thus the obtained relationships 
below hold when the corresponding maximum degree is max > 0. 

By the method of mathematical induction we can prove that: 


• For h( 2 a) • Maximum grade of polynomial K 2a (g) is max = q + a(m — 2) 
where a — 1,2,.... Maximum grade of the polynomial Z 2a (g) is 0. 


• For h( 2 a+i) ■' Maximum grade of polynomial of K 2a +i(g) is 0. Maximum 
grade of the polynomial of Z 2a+1 (g) is max = q + a{m — 2) — 1 where 
o'= 0,1,.... 
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4 


Examples 


4.1 Generalized Korteweg-deVries equation 

We shall consider the equation 


du du d 3 u 

m +Au + w 


(4.1) 


where A is a parameter and p is a positive integer number. This equation has 
different applications as for an example in the electrohydrodynamics [53]. For 
the case when p is a positive integer Eq. fld.lll is called generalized Korteweg- 

d 3 u 

deVries equation. It is obtained by addition of the dispersion term —— to 

ox 6 


the nonlinear convective wave equation ——K Au p — = 0. 

dt dx 

We search for solutions of Eq. (14.11) of the kind u = h[f(g(£))] where 
£ — ax + /3t, g is solution of the simplest equation (12.9|) and / is given by 
Eq. (l2.5|) . The substitution of u = h[/(g(£))] in Eq. (l4.1|) leads to equation of 
the kind (I2.2ip where 


W 0 (g) = uKi(g) + g,AK 0 (g) p K±(g) + g 3 K 3 (g) 

W\(g) = uZi(g) + /iAK 0 (g) p Zi(g) + g 3 Z^{g) (4.2) 


As K\ = Ji 3 = 0 there is no need to balance the relationship for Wo (g). The 
relationship for W\(g) has to be balanced however. The resulting balance 
equation is 

m = 2 + pq (4.3) 

Let us consider the case q — 1, m — 2 + p. Then from Eqs. (l2.5|) and (12.9[) 
one obtains 

2+p 

h = b 0 + big ; gf V) = ^ ajg J . (4.4) 

j =o 

Substituting corresponding relationships for the polynomials Z \, Z 3 and K 0 
in the second of equations (14. 2\\ we obtain the following system of nonlinear 
algebraic relationships among the parameters of Eq. d4.ip and the parameters 
of the solution: 

z'Mo, k+ (^) M^o~ fc& i +1 + ^p 3 (k + l)(k + 2)a k+2 bi = 0, k = 0,...,p. (4.5) 
where 6 is the delta-symbol of Kronecker. 
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Solution of Eqs. fl4.5p can be obtained when b 0 = 0. Then 

z'Mo.fc + +1 Sk )P + -p 3 (& + l)(fc + 2)ak+2b\ = 0, k — 0, 

The system (14.6[) is reduced to 


k — 0 : v + /i 3 a 2 = 0 

k = 1 : 03 = 0 


k = p — 1 : a p+ i = 0 

k — p : + ^p 2 (p + l)(p + 2 )a p+2 = 0 . 


The system (14. 71) has a solution: 


z/ 


a 2 — 


®p+2 


2A^ 


/i J ‘ p 2 (p+l)(p + 2 ) 

Hence the simplest equation Eg. (12.91) becomes 


9(i) 


v 


7 9 


2 b\ 


:9 


p+2 


pfiA p 2 A(p + 1 ) (p + 2 )' 

A solution of this simplest equation is as follows. The function 

9(f) = Jki 

(where co and fl are parameters) is solution of the equation 


Hence if 


then 


n 2 ,2 2 _ u 2+2 /u 

3(i) u 9 Q2/u}9 


2 2 p 2 AVt p l) p l 4p 3 

U p ^ 2 (p + l)(p + 2 ) ’ V p 2 


«(0 = 


Qb\ 


£ = px + ut 


cosh 2 /p ( 0 : 

is solution of Eg. (14.ID . 

Let us note that according to Eqs. fl4.10D and (14.111) 


K),0,4/p 2 ,0,...,a p+2 =--4^ 2 ,2 + P) - 


n 


cosh 2 / p (£) 


,p. (4.6) 

(4.7) 

(4.8) 

(4.9) 

(4.10) 

(4.11) 

(4.12) 

(4.13) 

(4.14) 
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Let us consider the particular case p = 1, A = —6. Then Eq. 04.ip is 
reduced to the classical Korteweg-deVries equation. Setting b\ — 1 and 
= —2 we obtain p — 1, v — —4. Thus the solution (14.13[) reduces to the 
one-soliton solution of the equation of Korteweg-de Vries 


u(x, t ) 


2 

cosh 2 (a; — At) 


(4.15) 


Let us stress the following. Above p was arbitrary non-negative integer. Now 
we shall show that p can be arbitrary non-zero real number. Namely we shall 
prove 


Proposition. Eq. 

value of p. 


H4.lt )|) is solution of Eq. fl 4-l\) for arbitrary real nonzero 


Proof. Let us substitute Eq. 04,1.31) in Eq. 04.ip where parameters a and /3 are 
given by Eq. 04.121) and p is an arbitrary nonzero real number. Eq. 04.11) is 
satisfied. Hence Eq. 04.131) is solution of Eci. 04.ip for arbitrary real p ^ 0. □ 

Let for an example p = 3/2. Then a solitary wave solution of the equation 


du 

dt 


+ An ^ + 

OX 


d 3 u 

dx 3 


0 


(4.16) 


is 


u(x,t ) = Qb\ cosh 4 / 3 



(4.17) 


4.2 Application of methodology to the higher order 
Korteweg-deVries equation 

Let us apply the above methodology to a more complicated equation such as 
the second-order Korteweg-deVries equation 


du du du du d 2 u 
dt dx + a ° u Q x + ai dx dx 2 


f a 2 u 


d 3 u 


,du 


\-0i3U~ — boq 
dx 6 dx 


d 3 u 

dx 3 


d 5 u 


0 (4.18) 


This equation is known also as Olver equation [53J and it is a second order 
equation for description of shallow water waves (the first order equation is 
the famous Korteweg-deVries equation). Eq. (14.18j) has been used also as 
a model equation for nonlinear waves in a liquid with gas bubbles [55l 56] 
From the point of view of the method of simplest equation (based on the 
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first step of a test for Painleve property) the equation was discussed in [55] - 
[57] and various solutions expressed by the elliptic function of Weierstrass 
have been obtained there. Below we shall discuss Eq. fl4.18fl from the point 
of view of the modified method of simplest equation (based on the concept 
for balance equation). We stress again that the modified method of simplest 
equation is a version of the method of simplest equation and then some of 
the obtained below solutions (especially for the case m — 3) will be the 
same as these obtained in m- After studying the methodologies from (S3 
and from this paper the reader will have an extensive understanding about 
the possibilities for obtaining exact analytical solutions of nonlinear partial 
differential equations on the basis of the method of simplest equation. 

We search for a solution of the kind u(x, t ) = w(£) = u(ax + /3t ) (We note 
that in [57] a = 1). According to the theory above Eq. (14.18H can be reduced 
to the following system of equations 

W 0 (g ) = 0 ; 

Wi(g) = (g + i')Zi(g) + a 0 gK 0 (g)Zi(g) + aig 3 K 2 (g)Z 1 (g) + 
a 2 g 3 K 0 (g)Z 3 (g ) + a 3 gK^(g) Z 3 (g) + a 4 g 3 Z 3 (g) + 
a 5 g 5 Z 5 (g) = 0 (4.19) 

The second of Eqs. fl4.19jl has to be balanced and the balance is as follows 

1. rn — 1: there is no balance; 

2. m — 2: there is no balance; 

3. m > 3: the equation is balanced if q = m — 2. 

Let us now consider several cases. 

4.3 Case m = 3 

In this case q — 1. Then 


u(0 = b o + M(0; i = gx + vt 

(4.20) 

9(i) — a o + a i 9 + a 29 2 + 0-39 3 

(4.21) 


In addition we have to solve the system of nonlinear algebraic equations 
for the parameters of the solution that can be obtained from the equation 
W\(g) = 0 from Eqs. fl4.19jl . 
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We note that the general solution of Eq. (l4.21l) is given by the special 
function 14 0jaija2jO3 (£; 1; 2, 3) and for the special case when a 2 = 0 and a 3 = 4 
we have reduction of V to the elliptic function of Weierstrass 


f // ao,Ql>0,4('C) 1) 2, 3) ®Cb ® 1 ) 


(4.22) 


Let us first consider the case of general non-reduced solution V a0!ai ,a 2 ,a 3 (£; 1, 2, 3) 
of Eq. (14.211) . The system of nonlinear algebraic equations arising from second 
of the Eqs. 04.191) is (IA.1D from the appendix A. The solution of the system 
dO) is: 


0-3 ~ 


O 2 — 


bi 


30/i 2 cc5 


^ — cxi — 2 a 2 + \J (oq + 2 a 2 ) 2 — 400 : 30 : 5 ^ 


5/i 2 a 5 + \J{oi\ + 2a 2 ) 2 - 40a 3 a 5 J - 

(a 2 b 0 + aq) y/ (aq + 2a 2 ) 2 - 40a 3 a 5 - aqa^o - 2«2&o + 
20a 3 cc5&o T 10ao a 5 — oqoq — c 2ot 4 oi2 


a 1 = — 


Ti 

T 2 


T, = - 


25/iala 3 a 5 bl — G/iay^o + 90fiala 3 a 5 bl + 25 / 10 : 00 : 40:560 + 

30 fia o ala 5 b 0 — 300/10:00:30:560 T 2 fiOi ^ oi 2 oi 4 bo — 2/10:40:20:460 T 

2 2 2 2 2 2 
50 / da o a 5 — 35^0:00:10:40:5 — 20 / ioiQd 2 d 4 Oi ^ T /loqoq — 6/10:20:4 T 

9 9 no o 

25/10:40:5 + 50/10:20:5 — 500 / iq : 3 q : 5 + 25110:40:5 + 50^0:20:5 — 

500 ziq : 3 q : 5 + 80/10:30:40:5 + 50/10:40:20:5 + 50 uaia 2 a 5 + fia 1 a 2 by — 

/laqa^o — 300fialalbl — l2fiala 4 bo — /ia: 10 : 20:4 — 20 / 10 : 40 : 20 : 30:560 + 

15/10:00:10:20560 — 70/10:10:30:40:560 T 120/10:20:30:40:560 T 

( 2 2 7 2 q,2 2 2 

6/iQ^2C^4&0 — /XQ^iQ^^o — 15 /XQ(oC^ 4C^5 “f“ — flCLiCi ^ -f- 3^Z&2^4 “1“ 

25 / 10 : 10:5 + 25 hq:iq :5 + 25 / 10 : 10 : 30:560 — 30 / 10 : 20 : 30:560 + 25 / 10 : 0010:560 — 

2/10:40:20:460 — 30/10:30:40:560 ) \ J (0:4 T 2 q: 2 )^ — 400:30:5 


T 2 = / 1 2 0:561 


(5o ? 2 T 3 a: — 2 2 — 300 : 30 : 5 )^/( 0:1 T 2q : 2 ) 2 — 40q: 3 q:5 T 


10q? 2 q:2 T 7o:4Q ? 2 — 1300 : 40 : 30:5 — 60:2 T 60d 2 d 3 O'^ 


(4.23) 
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The solution of the higher order Korteweg-deVries equation (j4.18j) is 

«(0 = bo + ftiVao.ai, 02 , 03 ^; 1,2,3); i = jix + vt (4.24) 


where ai i2j3 are given by Eqs. (14.2311 . 

Let us now consider the particular case (14.22j) when the V-function is 
reduced to the elliptic function of Weierstrass. In this case we have to set 
0,2 = 0; o 3 — 4 in the system (??). The solution of the obtained system of 
algebraic equations is 


bn = 


3/r 


b\ = 


OiQ\J (<3q Q! 2 )^ — 400305 — Q.q(X\ — 2ao«2 T 4oqCr 3 

2a 3 (oq + 2 qi 2 ) 2 — 40a 3 a 5 — oq^ 

^ — Oil — 2 oq + ■sj ( Oil T 2 ctq ) 2 — 40^305^ 


a 3 


(4.25) 


CL \ 

t 3 


T 3 

T,i 

1 

_ 3 


T 4 = 


0 0 0 o 0 0 

— — 2 /iC^Q(^x ^2 “ 1 “ 20 /iC^QQ^ 3 Q (5 — 8 /iQ(oC^ 2 ^ 3^4 “ 1 “ 

+4/oala 3 + 8/zaqa 2 a!3 + 8/ia\a 3 — SO/ia^aq + 4i/a^a3 + 8z/aqa 2 a!3 + 

8ua\a 3 — 80uala 3 + f/ialai — 4/iaqa 3 — 4z/aqaq'j \J (aq + 2aq) 2 — 40a 3 Q;5 

H 5 ^\/(oti + 2a 2 ) 2 — 40 o! 3 Q!5 - aq^ ^aq\/(aq + 2a 2 ) 2 - 40a 3 a 5 - 


— 2oqcr 2 V 12 ct 3 ct 5 ^ 

The solution of Eq. fl4.18p becomes 


OioW (oil + 2a 2 ) 2 — 40a 3 a5 — croQq — 2ao Q; 2 + 4aqaq 

m(0 =-7--- + 


3/r 


2a 3 ^\/ (oq + 2a 2 ) 2 — 40a 3 Oi 3 — oq^ 

( — oil — 2aq + \J (oq + 2a 2 ) 2 — 40cr 3 o;5^ 


a 3 


p(€',ao,ai) 


£ = fix + ut. 


(4.26) 
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where a± is given by the corresponding relationship from Eqs. fj4.25p . 

Let us now consider as simplest equation Eq. fl4.lip for the case cu — 2, 
namely 

g\ 1 ) = 4 g 2 - 4j s (4.27) 

In this case ao = cq = 0, 02 = 4 and a 3 = — (|. The solution of Eq. fl4.27p is 

9(f) = dbr wyeu. 2 .3) (4.28) 

One solution of the system of nonlinear algebraic equations for this case is 


&o —- - ^4 /j?O i\ T Afj ?012 + 

0 ) 3 y («i + 202) 2 — 400 ) 30)5 + 01 _ 

Oq) \J (03 T 202) 2 — 4 O 03 O 5 T 4p o^ T 12p 03 O 2 T 


8 yU^Q )2 — 80O)5yU 2 O3 T OqO} T 20q02 — 4 O 4 O 3 


(4.29) 


bi 

v 


3fi 2 

o 3 0 


^Oi + 2o 2 T \/(c>!i + 2 q?2) 2 — 400)305 


_h_ 

4o 3 (v^TT^F^lO^O! T o 2 T 203 O 2 T 2o^ — 2 O 03 O 5 ) 

^GqGi -b 4 q^ic^ 3 -b 16/i “I - 32/i g^G 2 “I - 16/i G1G2 — 

256/i 4 0i030s) \/(01 + 2 q ) 2 ) 2 — 400)305 + 64p 4 o'j ! a2 + 80/i 4 o 2 o 2 — 
576 /i 4 o 2 o 3 05 + 32/i 4 0i02 — 512/i 4 Q) 1 o 2 o 3 05 — 192fi 4 ala 3 a 5 + 1920 / u 4 a 2 o 2 — 

000 o 00 j o 

G 0 Gi — 2 GqG^G 2 “I - 2OGQG3G5 — 8G0G2G3G4 -f- SGgG^ + 4 g^G 3 ~b 8G1G2G3 -b 


80^03 — 800)305 


and the corresponding solitary-wave solution of Ea. fl4.18p is 

1 


u (0 ~ - h o NO An = - (4/i 2 Oi + 4/r 2 0 2 + 

o 3 V (“i + 2a 2 ) - 4O03O5 + 01 L 

Oo) y /(01 T 202) 2 — 4 O 03 O 5 T 4 n 04 T 12p O 1 O 2 T 

3 p 2 


8/Vo 2 — 80o 5 /i 2 o 3 T oqOi T 2 oq 02 — 40403 


+ 


03 


\/(01 + 2 o 2 ) 2 - 4 O 03 O 5 ) 


) cosh (0 


Oi + 2 o 2 + 
(4.30) 
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where 


4 « 3 (\/ («i + 2 a 2 ) 2 — 40a3«5ai + a 2 + 2 ai «2 + 2 ct 2 — 20 « 3 q; 5 ) 

/ n /IQ /10 /10 

(q^qC^i 4~ 4o^q?3 H - 16/x + 32/i 4~ 16/z GfiQ^ — 

256/i 4 aia30!5) \/(an + 2 q? 2) 2 — 40a3«5 + 64/i 4 afa2 + 80/i 4 a 2 a 2 — 

/IQ /IQ A A O /I o o 

576/i 0^0305 + 32fi aqa 2 — 512/r aqa^aqas — 192/i a 2 « 3«5 + 1920/i a 3 a 5 — 

OOO o 00 j o 

— 2c^QQfiQf2 ~h 200 ^ 0 ^ 3(^5 — 800 O 2 O 3 O 4 4“ 8 ogO^ 4~ 40^03 4~ 801 O 2 O 3 4~ 

8 a 2 a 3 — 80a 3 a 5 t 

4.4 Case m = 4 

In this case g = 2. Then 

w (0 = &o + 6 i 0 (O + 6 20 2 (O; ^nx + vt 

and 

fl'fi) = a o + «ig + a 2ir + «3g 4 + 1 

In addition we have to solve the system of nonlinear algebraic equations 
for the parameters of the solution that can be obtained from the equation 
Wi {(]) = 0 from Eqs. fl4.19p . 

The general solution of Eq.. (14.321) is given by the special function !4 0iaiia2ia3>a4 (£; 1, 2,4). 
For several particular cases this function can be reduced to the elliptic func¬ 
tions of Jacobi. Examples are 

Ei > 0 ,-(i + fc 2 ) A fc 2 (£;l, 2 , 4 ) = sn (£;&); 
hi-fc 2 ,o,( 2 fc 2 -i),o,-fc 2 (£, 1 , 2 ,4) cn(£, fc), 

^C(i-fe 2 ),o,( 2 -fc 2 ),o1; 2,4) = dn(£;fc). (4.33) 

For the general case the system of nonlinear algebraic equations becomes 
(IA.2D from the Appendix. This system of equations possesses exact solutions 
for 00 , 01 , 02 , 03,04 but they are very long and we shall not reproduce them 
here. Instead of this we shall give a solution for an illustrative particular 
case. From the first of equations (IA.2D one obtains a solution as follows 

(^2 [—c^i ~~ 2a 2 + \J (ot\ + 2a 2 ) 2 — 40a 3 a 5 ] 

4 120a,5/r 2 

Let us consider the particular case 

(o/.\ + 2 « 2) 2 

5 40ct 3 


(4.34) 

(4.35) 


(4.31) 

(4.32) 
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For this particular case a solution of the system of equations (IA.2|) is as 
follows 


o 4 

0 3 

0 2 

0\ 

Oo 


a 3 b 2 


3 /i 2 (of + 2 o 2 ) 
2 o 3 6 4 


3 /i 2 (o 4 -I - 2,012) 

4 :OiiOi 3 bob 2 4 ~ O4O364 4 - 200O462 4 ~ 400O262 — 8030462 
4 ^ 2620 : 1(01 + 2 o 2 ) 

61(1201036062 — 01O36 2 4- 6oqOi6 2 4 - I200O262 — 2403O462) 


12 /i 2 oi 62 (oi 4 - 02) 

336/iOiOg6o62 


84 /iOiOg 6 o 6 f 62 4 - 


48/i 3 of 0263(01 4 - 2 o 2 )( 7 oi — 6 o 2 ) 

7/1O4O364 — 288/10402036063 4 - 72/1O4O2O3606462 — 6/1O4O2O364 4 - 
336/10004036063 — 42/100O4O36462 4 - 384/100O4O2O36062 — 48/1O0O4O2O36462 — 
576/1O0O4O2O36062 4 - 72/100O4O2O36462 — I344/1O4O3O46062 4 - I68/1O4O3O46462 4 - 
II52/101O2O3O46062 — I44/101O2O3O46462 4 - 24/1O0O462 4 - 144 /io 0 O 4 O 26 2 4 - 
288/1O0O4O262 4 - I92/1O0O262 — 672/1O0O4O3O462 — I728/1O0O4O2O3O462 — 
768/100O2O3O462 4 - 2304/1O1O3O462 4 - 768/102O3O462 4 - 24O/1O4O362 4 - 


48O/10/O2O362 4 - 240^040363 4 - 48OHO/O2O362 


(4.36) 


and the solution of Eq. fj4.18j) for the discussed particular case f)4.35j) is 


u (0 —bo + biVa 0tai ,a 2 ,a 3 ,a 4 (^i 1,2,4) + b 2 V^ oaia2azaA (£; 1,2,4) (4.37) 

Let us now consider the particular case when V^ 0)ai)a2iO3iO4 (£; 1, 2,4) is 
reduced to the Jacobi elliptic function sn(£; k). In this case we have to set 
a 0 = 1; cq = 0; a 2 = — (14- k 2 ); a 3 = 0; a 4 = k 2 in Eqs. fi A.2 ll . The full solution 
of the system of algebraic equation is again very long. In order to illustrate it 
we discuss the particular case given by Eq. (I4.35|) . A solution of the obtained 
system is 


bo 


1 

2/10403(704 — 602) 


21 /i 3 o 3 4 - 24 /i 3 040 2 — 36 /i 3 0402 — 7 /j,a 0 al — 


8/100O4O2 4 - 12 /io 0 O2 4 - 2804/10304 — 24/io 2 o 3 o 4 4 - ^ — (704 4 - 8o 4 o 2 — 12 o|)/i x 
(21/i 5 04 4 - 24/i 5 o 3 o 2 — 36 /j, 5 alal — 5/1O0O4 4 - 20/iOqO2 — 80/ioo0 2 o 3 o 4 4 - 8O/1O3O4 4 - 
26/10403 4 - 2 OZ/O 4 O 3 ) 
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0 


h = 

b 2 =-—--- fi 3 ai(7al + 8aqa 2 — I2a^) + ( — (7a; 2 + 8a 4 a 2 — 12a 2 )/i x 

2paia 3 (7ai — 6a 2 ) |_ V 

(21yU 5 o4 + 24/i 5 a 3 a 2 — 3Qp^a\a\ — 3pa^a\ + 20/j.alal — 80/ia 0 a 2 a 3 a 4 + 80 pa\a\ + 
\ 1/2- 

20 pa\a 3 + 20z/a 2 a 3 ) 

k = —= -—--- ————-—-( ax(^ + 2a 2 )(7/n 3 af + 8fx 3 a^a 2 - 12fx 3 a!a 2 + 

\/ 2 /iQfx(c^i + 2 ^ 2 ) _/^(7o^i 00 ^ 2 ) \ 

(—(7a 2 + 8«ia 2 — 12a2)/i(21p 5 a4 + 24/i 5 a 3 a 2 — 36p 5 a 2 a 2 — 5/iaga 2 + 20 pa^a 2 — 

1/2 

(4.38) 


80/ia 0 a 2 a 3 a 4 + SO/xa^a 2 + 20pa 2 a 3 + 2Qva\a 3 )) 1 ^ 2 ) 


and the solution of Eq. (l4.18p is 
u(0 = 1 


21 /i 3 a 3 + 24/x 3 a 2 a 2 — 36/x 3 aia 2 — 7 pa 0 a 2 


2 noi\oi 3 (7 oi\ — 6a 2 ) 

8yuo;ottia 2 + 12paoa 2 + 28a 4 /ia 3 a 4 — 24pa 2 a 3 a 4 + f — (7a 2 + 8a 4 a 2 — 12 a 2 )/i x 

(21p 5 a 4 + 24p 5 a 3 a 2 — 36/x 5 a 2 a2 — 5/jalal + 20 pa 2 0 a 2 — 80pa 0 a 2 a 3 a 4 + 80 pa 3 a 2 + 

3 


20pa 2 a 3 + 20z/a 2 a 3 ) 


1/2 


/i 3 a 4 (7a 4 4~ 8a 4 a 2 — I 202 ) + 


2/zaia 3 (7ai — 6a 2 ) 

(7a 2 + 8a 4 a 2 — 12a 2 )p x (21/i 5 a 4 + 24/i 5 a 3 a 2 — 36/i 5 a 2 a2 — 5/zaga 2 + 


9 9 99 9 9\\ 

20/ia 0 a 2 — 80/xaoa 2 a 3 a 4 + 80/xa 3 a 4 + 20 / ua 1 a 3 + 20z/a 4 a 3 ) ) 


\ 1 / 2 - 




X 


sn 2 s f; - 7 =- 7 --—- —pz - - —r ( ai(«i + 2a 2 ){7p 3 a\ + 8 /i 3 a 2 a 2 - 

V V2/iai(cti + 2a 2 ) Lh(7«i - 6a 2 ) V 

12yU 3 aia| + (—(7a 2 + 8 a 4 a 2 — 12 a 2 )/i( 21 p 5 a 4 + 24/x 5 a 3 a 2 — 36p 5 a 2 a 2 — 

0 0 0 0 0 0 0 9\\1 / 0 \ 

5/xaga 4 4~ 20 /ictgCt 2 — 80paoa 2 a 3 a 4 4“ 80/xa 3 a 4 4- 20pa 4 a 3 4~ 20Z 2 a 4 a 3 )) t 
£ = px + id. 

Let us now consider the particular case when the simplest equation is the 
Riccati equation 


1/2 


(4.39) 


<7(i) — c o + c i 9 + c 2 g' 


In this case 


®o — Cq; cl\ — 2coc 4 ; ci 2 — c 2 + 2cqc 2 ; a 3 — 2c 4 c 2 ; a 4 — c 2 


(4.40) 

(4.41) 
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The full solution of the system of algebraic equation is again very long. We 
illustrate it for the particular case given by Eq. (l4.35p . A solution of the 
system of algebraic equations is 


bo 

bi 

b 2 


8 /i 2 a^coC 2 + n 2 a\c\ + 16 p 2 oi 02 CoC 2 + 2 p 2 aq a 2 c 2 + 2 ooOi + 4 a 0 «2 — 80304 

4oio 3 

3/i 2 CiC2(oq + 202 ) 

ot 3 

3p 2 c 2 (aq + 2a 2 ) 
a 3 

— — ( 12 p 4 a^CoC 2 — 56/i 4 a 4 c 0 CiC 2 + 7 p 4 o 4 c 4 + 128/i 4 afa 2 CoC2 — 

80afa3 V 

4 q n /i q /| A O O O O A O O O 

64 fi a 1 a 2 cocJC 2 + 8/i a 1 « 2 C 1 — 192p a 1 a 2 C 0 c 2 + 96p a 1 a 2 c o c i c 2 — 
12 p 4 a 2 a 2 C 4 — 20«otti + 80 oqO2 — 320 ooo 2 o 3 04 + 320ogO 2 + 80o 2 o 3 ] (4.42) 


From the solutions of the Riccati equation (14.391) here we shall discuss the 
solution (c is a constant of integration) 


s (0 = 



(c? 


4c 0 c 2 ) 1/2 

2 


(Z + c) 



(4.43) 


valid when c 2 > 4coC2 and ( 2c2g (£) + C 1) 2 < c 2 — 4coC2- Then the solution of 
Eq. (14.181) is 


«(0 = 


8p^a 2 coC2 + /tro 2 cf + 16 /x 2 oi 02 CoC 2 + 2)Lroi02Cj + 2a 0 «i + 4 oq 02 — 80304 


tanh 


3p 2 ci(oi + 2 o 2 ) 

2o 3 

_ 3^t 2 («i + 2 q 2 ) J tmh 
4 o 3 ( 


4 oio 3 

(ci - 4cqc 2 ) 1/2 


(Z + c) 


(4 - 4c 0 c 2 ) 1/2 


(£ + c) 


~ ci r 

2 


- Cl 


(4.44) 


which is a kink. 


4.5 Case m = 5 

In this case q = 3. Then 

u{Z) = + &2</ 2 (0 + 63^ 3 (05 f W + ut 

and 

,9p) — a o + a i9 + « 2 lT + a 3 9 3 + ° 49 4 + a 59 5 


(4.45) 

(4.46) 
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In addition we have to solve the system of nonlinear algebraic equations 
for the parameters of the solution that can be obtained from the equation 
W\(g) = 0 from Eqs. fl4.19j) . 

The general solution of Eq. (I4.46j) is given by the function Va 0 , ai ,a 2 ,a 3 ,a 4 ,a 5 (^', 1, 2, 5). 

The system of algebraic equations obtained from Eq. fl4.18p is (1A.3I) . The 
full solution of this system is quite long. We shall illustrate the solution for 
the particular case 

(o 4 + 202) 2 (oq + 2) 2 


«2 — «3 — 1 ; Ot! 5 — 


40aq 


40 


(4.47) 


The solution for this case is 
46 3 


a 5 — 


<24 = 


<23 — 


o 2 — 


81/i 2 

206 2 

243/x 2 

4(6 2 + 76 4 6 3 ) 

243/i 2 6 3 

2(162oi6 0 6 2 + 72 oi6i6 2 6 3 — 10aq&2 + 243a 0 ^3 — 324o 4 6 2 ) 

2187/i 2 6 2 Oi 


a 1 = 


T 


6561 /i 2 aq &3 \ 


( 162oi 6 0 6 2 6 2 + 108aq& 2 &3 — 6301616263 + 8aq&2 + 


243«o^2^3 


324o46 2 6: 


ao — 


6561 /i 2 aq &3 


648ai&o&i&3 — 162ai6 0 &2^3 — 144 oi6 2 626 2 + 6801646263 — 80163 + 


972006163 — 243006363 — 1296046163 + 324o 4 6 2 6 2 


/2 = 


— 656 IOZ/O 163 

. 2r 2;,4 


= 4374of6g63 — 2916 o 2 6 0 6i6 2 63 + 64801606263 — 216 o 2 6 f &3 + 702 o 2 6 2 6 2 6 2 - 
288o 2 6i6'26 3 + 32 o 2 6 ® + 13122o 0 oi6 0 63 — 43 74o 0 oi6i6 2 63 + 97200O1626 2 
17496 oiO 4 6 0 63 + 583201O4616263 — I29601O4626 2 + 59049 oo &3 — 
288684o 0 o 4 6^ + 279936 o^ + 65610 oi 6 * 

and the solution of Eq. fl4.18j) (for the particular case given by Eqs. fl4.47p ) is 


(4.48) 


u{0 = 6 0 + 61 K 0 , 


ai, 0 , 2 ,< 13 , 014 ,as 


(£f 1 ; 2 , 5 ) + 6 2 V a 


ao,ai ,a2,a3,a4,a5 


63 K 


ao,ai,a2,<i3, a 4,“5 


(f; 1 , 2 , 5 ) 


(f; i,2,5) + 

(4.49) 
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where 


. 65610z/ai6o 

£ =-^——£ + vt 

^ 'T 1 * 


(4.50) 


4.6 Case m = 6 

In this case g = 4. We note that with increasing m (and q) the number of 
equations in the nonlinear algebraic system we have to solve becomes large 
very fast. When the number of equations become larger that the number 
of parameters of the solution then the system could not have any nontrivial 
solution. 

For the case m = 6 

w(0 = bo + big{£) +b 2 g 2 (£) + hg 3 (^) + 6 4 0 4 (£); £ = vx + vt (4.51) 

and 

9(i) — a o + Oi 9 + a 2 9 2 + g 3 + g 4 + OSS ' 5 + (4-52) 

In addition we have to solve the system of nonlinear algebraic equations 
for the parameters of the solution that can be obtained from the equation 
W\{g) = 0 from Eqs. fl4.19j) . The general solution of Eq. (14.521) is given by the 

function V^Q ,01,02,03,04,a.5 ,ae('C) 1 ) 2 , 6 ). 

The system of nonlinear algebraic equations connected to Eq. fl4.18|) is 
very large. It still has a solution. In order to illustrate this solution we shall 
consider a particular case of Eq. fj4.52p namely the equation of Abel of first 
kind 

0(i) = c 0 + ci <7 + c 2 g 2 + c 3 g 3 (4.53) 

The corresponding system of nonlinear algebraic equations consists of 12 
equations. The solution is very long and we shall illustrate it for the partic¬ 
ular case when a 5 = (ai 4 t 2ct2 . For this case one solution of the system of 
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algebraic equations is 

4/r 2 C2(ai + 2 ck 2 ) 


b A = 
bs = 
b 2 = 
bi = 
bo = 


3 a 3 cl 

16/^ 2 c|(ci!i + 2a 2 
3q<3 C\ 

S/i 2 c\{ai + 2a 2 ) 
cn 3 


16 / u 2 c 2 (3aic 0 c 2 + 4«ic 2 + 6 a 2 c 0 c 2 + 8 a 2 c 2 ) 

I 503 C 3 

00 000 o 00 

96/i aqc 0 c 2 + 8/1 aqcq + 192/i aqa^co^ + IQjji a\a 2 c 1 + 


c 3 = 

c 2 = 

» = 


30«ia 3 
15ao«i + 30ao« 2 — 60ct 3 a4 


27c 2 


3Cn 


4z/a 2 a 3 


a^af — 4aga| + 16aoa 2 asa4 — 16a|a| — 4a 2 «3 
The Abel equations becomes: 


(4.54) 


Jd) - C„ + c ig + J_J 2 + 27cl 


‘ 9 3 


Its solution is: 


9(0 = V 


5 cf 20 c? 5 c? 2 c? 


9 _ oc, ^uc 1 

Cn,2coci,-=l loyAj 


(4.55) 

(4.56) 


The solution of Eq. (14.18|1 is 


w(0 = & o + & 1 E 


c? 20c? 5 c? 2c? c? + &2E 


c 2 2cnCl m Jli U , ■ “ ' 2 2 rnr . Ai Ati 5c i 2c i c i , 

0 ’ 0 1 ’ 3 ’ 27 c 0 ’ 27 cq ’ 81 cg ’ 729 cg C 0 , 2 c 0 Cl, 3 , 27 c 0 ’ 


+ 


b , V 3 


5 c? 20c? 5 c? 2c? c? T b^V gc 2 2 Q C ? 5 c? 2c? 


r>2 O^^n III 1 1 1 II n 2 „ III II 

O’ 2 0 !’ 3 ’27c 0 ’i^3-’iI3-'7i5^ C 0 ,2c 0 Ci, 3 , 27 ,„, 


(4.57) 


where 


f = 


4i/a 2 a 3 


3 ’ 27 c 0 ’ 27 cq ’ 81 c 3 ’ 729 cg 


X + ut 


CKgCK 2 — 4aQa 2 + 16aott2Q; 3 a4 — 16a|a| — 4afo:3 
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2 ci 


9c3 


Let now Cq — ( C \ — ^ ) and C3 < 0. Then the Abel equation 


C 2 ( 2 Co 2 o 

9(1) = 3^ l' Ci ~ 5^ + Ci9 + 929 + 039 


(4.58) 


has the solution [58 


9(0 = 


exp 


* - 4 ) * 


C — 2 c 3 exp 




C 2 

3 c 3 


(4.59) 


where C is the constant of integration. Let us assume that C > 0. The 
use of the Abel equation f)4.58j) leads to a following solution of the system of 
nonlinear algebraic equations (for the particular case a 5 = ( “ 1 4 ^ 2 ' > ) 


b d = 


b-i = 


62 — 


61 = 


v = 


12 c^(oi\ T 2 ct 2 ) 
ct 3 

16/i 2 c 2 c 3 («i + 2« 2 ) 
a 3 

4/r 2 (3aiCiC 3 + c 2 + 6a 2 CiC 3 + 2 a 2 c|) 
a 3 

8yU 2 (9aqcic 3 — aqc| + 18a 2 cic 3 — 2 a 2 c|) 


/i 


9a 3 c 3 

2268^ 4 aqc 4 ^ — 3024/r 4 aqc 3 c 2 C3 + 1512yU 4 a 4 c 2 C2C 2 — 


1620afa 3 c| _ 

336/i 4 a 4 ciC2C 3 + 28// 4 aqC2 + 2592/i 4 a 3 a 2 c 4 c 4 — 3456yU 4 a 3 a 2 c 3 C2C3 + 
1728/i 4 a[l] 3 a 2 c 2 C2C 2 — 384/i 4 a 3 a 2 CiC2C3 + 32/T 4 — 3888/i 4 a 2 a2C 4 Cg + 

5184/i 4 a 2 «2C 3 C2C3 — 2592yU 4 a 2 a2C 2 C2C 2 + 576/i 4 a 2 a 2 CiC2C 3 — 48yU 4 a 2 a|c2 — 

405 ao a i c 3 + 1620 aQa 2c 3 — 6480 a 0 a 2 « 3 « 4 C 3 + 6480 « 3 a 2 C 3 + 1620 a 2 a 3 C 3 

(4.60) 
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Hence the solution of Eg. (14 .18 II becomes 


, 8/i 2 (9aiCiC 3 - «ic| + 18a 2 cic 3 - 2a 2 c|) 
«(£) = »o---x 


( 


exp 


9a 3 c 3 

c i~CK 


\ 


c 2 


yjc - 2c 3 exp 
4/i 2 (3aiCiC 3 + c 

2 ( c ‘ - il) « 

*iC2 + 6a 2 cic 3 
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(4.61) 


that is a solution of kink’s kind as the solution il 1. I III . 


5 Concluding remarks 

In this article we have discussed a version of the method of simplest equation 
applicable to a class of nonlinear partial differential equations that are much 
used as model equations in the area of natural sciences. Eq. (12.9)1 was used as 
simplest equation and we have described a methodology based on the concept 
for balance equations. This methodology reduces the studied nonlinear par¬ 
tial differential equations to systems of nonlinear algebraic equations. Any 
nontrivial solution of the obtained system of algebraic equations leads to an 
exact solution of the corresponding nonlinear partial differential equation. 
Discussed examples have shown the effectivity of the methodology of this 
variant of the method of simplest equation. 
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A Systems of algebraic equations connected 
to Eq. ( 14.181 ) 

A.l Case m = 3 

0 — — 04/1^0364 4 " 4 ~ O36 2 4 " -^“Os/^O 2 

0 — 15 /i 4 020 3 05 4 “ /i 2 o 2 Oi64 4 ~ /A02O264 4 ~ 3 /i 2 o 3 026 o 4 ~ 8/1^0304 4 - 
2036064 4 - 0064 

0 = /i + /i + o 0 /i6 0 + -ai/( 3 ai^i 4- a2fJ 3 b 0 a 2 + o 3 /i6q + 

o 4 /i 3 a 2 + a^/j 5 ^a 3 ai 4 - a]A (A.l) 

A.2 Case m = 4 


0 = 360 /i 4 o 2 o 5 4 - 6/i 2 a 4 Oi62 4 - 12/i 2 o 4 o 2 62 4 - o 3 6 2 
0 = 2 ai/i 2 b- 2 (50362 + 20461) 4 - 6oi/i 2 o 4 &2&i + 2Aa 2 /j 2 bia4 : b 2 + 

02/i 2 6 2 (160362 A 60464) 4 ~ 5 o 3 6 i 6 2 4 ~ o 5 /i 4 (840a3a 4 62 4 ~ 120 o 2 6 i) 

0 = 2 o 0 6 2 A 201/^(^(40262 + (3/2)0364)62 4 ~ 04/1^(50362 4 ~ 20464)64 4 - 

24o2/i 2 6 0 a 4 62 + 02/i 2 6i(15a 3 62 + 60464) + a 2 ^i 2 b 2 (8a 2 b 2 + 3 a 3 6 i) + 
203(26062 4 “ 6 2 )6 2 4 ~ 2 o 3 6 2 6 2 4 ~ 2404/(^0462 4 * O5/1 4 (48O02O462 4 “ 2IO0362 4 ~ 
120030464) 

0 = 3 oo 6 i 6 2 4~ 2 oi / i 2 (3 o 46 2 4- 0264)62 4- oi / i 2 (4 o 262 4- (3/2)0364)64 4~ 

02/1^69(150362 4 ~ 60464) 4 - 02/1^64(80262 4 ~ 30364) 4 “ cx 2 /j, 2 b 2 (80462 4 ~ 0264) 4 - 
4 o 3 6 0 6 i 6 2 4 - 0364(26062 4 - 6 2 ) 4 - o 4 /( 2 ( 15 a 3 6 2 + 60464) 4 - o 5 /i 4 [( 45 / 2 )a 2 6 i 4 - 
270 a 4 62 Oi 4 - 195036202 4 - 60046402] 

0 = 2 (zi + /1)62 4 ~ 2 o 0 /( 6 0 6 2 4 ~ O0/16 2 4 - 2oi/i 3 (2oo62 4 - (1/2)0464)62 4 - 

oi/i 3 ( 3 ai 6 2 4 - 0264)64 4 - a 2 /i 3 b 0 (8a 2 b 2 4 - 3 o 3 6 i) 4 - o 2 /i 3 6 i( 3 oi 6 2 4 - a 2 6i) 4 - 
2 o 3 /i 6 q 6 2 4 " 2 o 3 /i 6 0 6 2 4 - O4/1 3 (8a, 2 b 2 4 - 30364) 4 " 05/i°(144oo046 2 4 ~ 

84040362 4 ~ 3O04O464 4 ~ 320^62 I502O364) 

0 = (zi 4 " /i)6i 4 " 00/16064 4 “ Oi/i 3 (2ao62 + (1/2)0464)64 4 ~ o 2 /i 3 6o(3ai62 + o 2 6i) 4 ~ 
O3/16064 4 ~ 04 /i 3 ( 3 oi 62 4 - o 2 6i) 4 ~ Os/i^g^i 4 ~ I20464O0 4 ~ I50262O4 4 ~ 
3O0362O0 4 ~ (9/2)036404) (A.2) 


27 




A.3 Case m = 5 


0 = (81/2)aifi 3 a 5 bl + 81a2^ 3 bla 5 + 3a 3 nbl + (10935/2)a 5 ii 5 alb 3 
0 = oi/i 3 ( 12 o 4 6 3 + 7a 5 b 2 )b 3 + 27a 1 fi 3 a 5 b 3 b 2 + 81a 2 fi 3 b 2 a 5 b 3 + 

a 2 li :i b 3 (6Qa 4 b:s + 35a 3 b 2 ) + 803/16263 + a 3 /i 5 (7656a 4 a 3 b 3 + 1820alb- 2 ) 

0 = 3a 4 /i 3 ((21 / 2) a 3 b 3 + 6o 4 6 2 + (5/2)0563)63 + 203 / 1 3 (120463 + 7ci 3 b 2 )b 2 + 
(27/2)a 1 /i 3 a 5 b 3 bi + 81o 2 /i 3 6ia56 3 + o 2 /i 3 62(60a 4 63 + 35 a 5 6 2 ) + 

Oi 2 /-i 3 b 3 (420363 + 24 o 4 &2 + lOa.563) + 30:3/1(26363 + 62)63 + 403/16263 + 

03/16363 + O5/1 5 (25200463 + 3850563 + (10605/2)056303 + 2394056204) 

0 — 300/163 4 - 3o3/i 3 (9a 2 63 4 ~ 5 o 3 6 2 4 - 20463)63 4 - 2oi/i 3 ((21/2)o363 4 ~ 60462 
(5/2)0561)62 + oi/i 3 (12a 4 6 3 + 70562)63 + 81 o 2 /i 3 6 0 O563 + 
o 2 /i 3 6i (600463 4 ~ 350562) Qt 2 /i 3 b 2 (420363 4 “ 24o 4 62 4 ~ IO0563) 4 ~ 

02 /i 3 6 3 (27 o26 3 4 - 15 a 3 6 2 4 - 60463) 4 - 3o 3 /i(26 0 6 3 4 - 26 i 6 2 ) 6 3 4 - 2 o 3 /i( 26 i 6 3 4 - 62)62 + 
203/1626363 4 ~ 81o4/i 3 Q'563 4 ~ 05/1^(3645020563 4 ~ 324O03O463 4 ~ 

I56O03O562 4 ~ 72O0462 4 " 46204O563) 

0 = 5 o 0 /i 6 2 6 3 + 3 oi/i 3 (( 15 / 2 )ai 6 3 + 4 a 2 6 2 + (3/2)0363)63 + 2 oi/i 3 ( 9 a 2 6 3 + 

50362 4 " 20463)62 + 03/i 3 ((21/2)a 3 63 4 ~ 60462 + (5/2)0563)63 + 
o 2 /i 3 6 0 (600463 4 ~ 350562) 4 ~ o 2 /i 3 63(42a 3 63 4 ~ 24 o 4 6 2 4 ~ IO0563) 4 ~ 
a 2 /d 3 b 2 (27 a 2 b 3 + 15 a 3 6 2 4 - 60463) 4 - o 2 /i 3 6 3 (15ai63 4 - 8a 2 6 2 + 3 a 3 6 i) 4 - 
303/1(26962 4 ~ 64)63 4 - 203/1(26963 4 - 26362)62 + 03/1(26363 + 62)63 4 - 
04/1 3 (600463 4 - 350562) 4 - 05/1 5 (2490056303 4 - 2040046302 4 -1015056202 4 - 
9450363 4 - 840046203 4 - 1200463 4 - (555/2)056303) 

0 — 409/16363 4 ~ 200/169 4 ~ 3o3/i 3 (6a 0 6 3 4 ~ 3 oi 6 2 4 ~ 0263)63 4 ~ 2o3/i 3 ((15/2)o36 3 4 " 

4 o 2 6 2 4 - (3/2)0363)62 + oi/i 3 ( 9 a 2 6 3 + 5a 3 b 2 4 ~ 20463)63 4 ~ o 2 /i 3 6o(42o363 4 ~ 

24 a 4 6 2 4 - 10 a 5 6 i) 4 - o 2i /1 3 63( 27 a 2 b 3 + lba 3 b 2 + 60463) 4 - 02 /i 3 6 2 ( 15 ai 6 3 4 - 8a 2 6 2 4 - 
3 a 3 6 i) 4 “ a 2 n 3 b 3 ( 6 a 0 b 3 4 ~ 3 ci\b 2 + o 2 6i) 4 ~ 603/1696363 4 ~ 203/1(26962 4 ~ 63)62 4 ~ 
03/1(26963 4 - 26362)63 4 - 04/i 3 (42a 3 63 4 - 24 a 4 6 2 4 - 100563) 4 - O5/i 5 (1680a,00563 4 - 
I26OQ3O463 4 ~ 66O03O562 4 - IO5O02O363 4 - 48OQ2O462 4 ~ I7O02O563 4 " 2IO0362 4 ~ 
120030463) 



— (3 (y + /i))&3 + 3 aofib 0 b 3 + 30:0/16462 4 ~ 3 oi \ p ?( 2 aob 2 4 - (1/2)0161)63 + 

2o4/x 3 (6a 0 6 3 + 3 ai 6 2 + ^261)62 + oi/i 3 (( 15 / 2 )ai 6 3 + 4 a 2 6 2 + 

(3/2)0361)61 + a 2 ^ 3 b 0 ( 27a 2 b 3 + 150362 + 60-461) + 02/ i 3 6 i (150463 + 80262 + 3ci 3 b \) + 
oi2jJ?b2 (60-063 + 3 a-i 6 2 + 0-264) + 30:3/16063 + 403/1696162 + 03/1(26962 + 61)61 4 - 
o 4 /i 3 (27a 2 63 + 150362 + 60461) + 05/1 5 (243a|63 + (45/2)0361 + 756046309 + 

42O0562O0 4 ~ (1107/2)03630-4 + 27O0462O1 4 “ IO50561Q4 H - I950362O2 4 - 60046102) 

0 = ( 2 (z/ + /i))6 2 + 209/16962 4 ~ 09/163 4 ~ 2oi/i 3 (2a- 0 62 + ( 1 / 2 )oi 6 i) 6 2 4 - 
oi/i 3 (6a 0 6 3 4 - 3 a-i 6 2 + 0264)64 + O 2 /i 3 6 0 ( 15 a-i 6 3 4 - 8a 2 6 2 4 - 3 a 3 6 i) 4 - 
02/i 3 6i (60-063 4 ~ 3 a-i 6 2 + a 2 6i) 4 - 203/16062 + 203/16963 + 04/i 3 (15oi63 4 - 80262 4 ~ 3a 3 b \) 4 * 
05/1 5 (270090363 4 - 144090-462 4 - 60090561 4 -1950-10263 4 - 840-10362 4 - 30010461 4 - 
320-262 4 - 15020361) 

0 = {y 4 ~ /i)6i 4 ~ 09/16961 + Oi/i 3 ( 2 a 0 6 2 + (1/2)0-464)64 4 - o 2 /i 3 69(60-963 + 3 a-i 6 2 4 - 0264) 4 " 
03/16961 4 - 04/i 3 (6o963 4 - 3 a-i 6 2 + ( i 2 bi ) 4 - 05/i 5 ((45/2)a-l63 4 - 0-261 4 " (9/2)03610-1 4 - 
1502620-1 4 - 12046109 4 - 600-26309 4 - 30036209) 

(A. 3 ) 
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